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Abstract 

Recent financial disasters have emphasised the need to accurately predict extreme financial losses 
and their consequences for the institutions belonging to a given financial market. The ability of 
econometric models to predict extreme events strongly relies on their flexibility to account for 
the highly nonlinear and asymmetric dependence observed in financial returns. We develop a 
new class of flexible Copula models where the evolution of the dependence parameters follow a 
Markov-Switching Generalised Autoregressive Score (SGASC) dynamics. Maximum Likelihood 
estimation is consistently performed using the Inference Functions for Margins (IFM) approach 
and the Expectation-Maximisation (EM) algorithm specifically tailored to this class of models. 
The SGASC models are then used to estimate the Conditional Value-at-Risk (CoVaR) and the 
Conditional Expected Shortfall (CoES), measuring the impact on an institution of extreme events 
affecting another institution or the market. The empirical investigation, conducted on a panel 
of European regional portfolios, reveals that the proposed SGASC model is able to explain and 
predict the evolution of the systemic risk contributions over the period 1999-2015. 

Keywords: Markov-Switching, Generalised Autoregressive Score, Dynamic Conditional Score, 
Risk measures, Conditional Value-at-Risk, Conditional Expected Shortfall. 


1. Introduction 

Recent financial disasters have emphasised the need to accurately predict extreme financial losses 
and their consequences for the institutions’ financial health and, more generally, for the safety of 
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the broader economy. Major financial crisis, such as the Global Financial Crisis (GFC) of 2007- 
2008 and the European Sovereign Debt Crisis (ESDC) of 2010-2011, usually spread over the whole 
economy leading to sharp economic downturns and recessions. In fact, during huge crisis episodes, 
the failure of banks and financial institutions is not rare and may trigger other non-financial insti¬ 
tutions through the balance sheet and liquidity channels, threatening the stability of real economy, 


see, e.g., Adrian and Brunnermeier (2014), Adrian and Shin (2010), Brunnermeier and Pedersen 


(2009), [Brunnermeier et al. (2009) and Brunnermeier (2009). The ability of econometric models 
to predict such extreme events strongly relies on their flexibility to model the highly nonlinear and 


asymmetric dependence structure of financial returns, see e.g., McNeil et al. (2015). Despite its 
widespread use, the simple linear correlation fails to capture the important tails behaviour of the 


joint probability distribution, see, e.g., McNeil et al. (2015) and Embrechts et al. (1999 2002). 
Hence, modelling the tail dependence and the asymmetric dependence between pairs of assets have 
been becoming increasingly more important in a multivariate environment, especially after the 
recent crisis episodes. The departure from the linear correlation as measure of dependence usu¬ 
ally implies to go beyond the multivariate Elliptical assumption for the joint distribution of asset 
returns. In this respect, the copula methodology allows to model a huge variety of dependence 


structures, see, e.g., Durante and Sernpi (2015). Another interesting feature of the dependence be¬ 
haviour of stock returns, which displays its relevant role when measuring extreme co-movements, 
is that it usually evolves smoothly over time as a function of past assets co-movements, see, e.g., 


Engle (2002) and Tse and Tsui (2002). Due to the exposure of common shocks affecting all the 


market participants, the conditional correlation between asset returns increases during periods of 


financial instability, see, e.g., Kotkatvuori-Ornberg et al. (2013), Sandoval Junior and De Paula 


Franca (2012), Syllignakis and Kouretas (2011) and Kenourgios et al. (2011). Dynamic copula 


models are referred to Patton (2006) and Jondeau and Rockinger (2006), although the problem of 


modelling the joint co-movement of stock returns was already present in Bollerslev et al. (1988), 


Bollerslev 


and Engle et al. (1990), among others. Moreover, occasionally, we observe breaks 


into the dependence structure, which are more evident during crises periods and other infrequent 
events, as documented, for example, by Bernardi et al.| (2013b) and Bernardi and Petrella (2015). 
As regards dependence breaks, Markov switching (MS) models have been proven to effectively 


capture non-smooth evolutions of the volatility and correlations dynamics. Chollete et al. (2009) 
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and Rodriguez (2007), for example, have firstly adopted MS copula with static regime-dependent 
parameters to analyse financial contagion. 

In this paper, we propose to model the regime dependent dynamic of the copula parameters 


using the score driven framework recently introduced by Creal et al. (2013) and Harvey (2013). 
Specifically, we allow the copula dependence parameters to depend on the realisations of a first 
order Markovian process with a specific Generalised Autoregressive Score (GAS) dynamic in each 
regime, while retaining an appropriate arbitrary specification for the marginals’ conditional dis¬ 
tribution dynamics. In this way, we extend the GAS literature by introducing nonlinearity and a 


stochastic evolution of the dependence structure similarly to Boudt et al. (2012). We name this 
new class of models Switching Generalised Autoregressive Score Copula (SGASC) models. During 
last few years, the conditional score approach has been extensively used to model the time varying 


behaviour of unobservable parameters in an observation driven environment (Cox et al., 1981). The 
use of the score as general updating mechanism has been justified in several ways by the literature. 


Harvey (2013), for example, represents the score process as a filter of an unobservable component 


model, while Creal et al. (2013) suggest that the use of the score to update the latent parameter 
dynamic, can be interpret as a steepest ascent algorithm for improving the model’s local fit given 
the current parameter position, as it usually happens into a Newton-Raphson algorithm. Recently, 


Blasques et al. (2015) and Blasques et al. (2014b) have demonstrated that the score driven processes 


are optimal within the class of nonlinear autoregressive dynamic. More precisely, they argue that, 
only the GAS processes are optimal in the sense of reducing the local Kullback-Leibler divergence 
between the true and the model implied conditional densities. Furthermore, the Kullback-Leibler 
optimality property of GAS models holds true under very mild conditions irrespective to the level 
of the possible model misspecification. Several theoretical results for the maximum likelihood es¬ 


timate of score process has been developed by Blasques et al. (2014c), Andres (2014), Blasques 


et al. (2014a) and Harvey (2013). Moreover, score driven processes have been proved to be ef¬ 


fectively used in many empirical applications. Most applications cover volatility modelling, as for 


example, Harvey and Luati (2014), Harvey and Sucarrat (2014), Caivano and Harvey (2014) and 


Creal et al. (2011a). Other empirical applications are in the systemic risk measurement as, for 


example, Blasques et al. (2014d), Lucas et al. (2014a) and Oh and Patton (2013), in the credit risk 


analysis, Creal et al. (2011b), in macroeconomics, Massacci (2014) and Bazzi et al. (2014), and in 
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dependence modelling, Harvey and Thiele (2014), Janus et al. (2014), and De Lira Salvatierra and 


Patton (2015). 


The documented superior ability of GAS filters to approximate complicated nonlinear data 
generating processes in a straightforward and effective way (see, e.g., Koopman et al.][2015) is par¬ 
ticularly useful in the context of dependence modelling through copulas here considered. Specif¬ 
ically, the use of score driven models really helps to deal with cases where it is not clear how to 
update the parameter dynamics as for the archemedian copulas. Another relevant contribution of 
this paper is to introduce and estimate the Conditional Value-at-Risk (CoVaR) and the Condi¬ 
tional Expected Shortfall (CoES) risk measures, recently proposed by Adrian and Brunnermeier 


(2011, 2014) and Girardi and Ergiin (2013), for the class of SGASC models. The CoVaR measures 
co-movements between any two distinct institutions by extending the Value-at-Risk (VaR) to a 
conditional approach. Following the CoVaR methodology, the risk of an institution is evaluated as 
its VaR conditional to a relevant extreme event affecting another institution. The appealing charac¬ 
teristic of the CoVaR risk measure is that it inherits the flexibility of the dynamic switching copula 
framework here developed. The copula approach naturally adapts to environments characterised 
by different kind of upper and lower tail dependence enabling the CoVaR as an effective measure of 
the extreme conditional co-movements among financial variables. The literature on co-movement 


risk measures has proliferated during the last few years, see, e.g., Bernardi et al. (2015), Bernal 


et al. (2014), Castro and Ferrari (2014), Girardi and Ergiin (2013), Jager-Ambrozewicz (2013), 


Sordo et al. (2015). Bisias et al. (2012) provide an extensive and up to date survey of the systemic 


risk measures that have been recently proposed. 

One of the main appealing characteristics of the copula framework, with respect to standard 
distributions, relies on its ability to model the marginals’ dynamics separately from the joint de¬ 


pendence structure, see, e.g., Nelsen (2007). Marginals and dependence separability has some 
additional advantages even from the econometric point of view, since it permits to employ a two- 
step procedure to estimate the parameters. This two-step procedure is known as Inference Function 


for margins (IFM) and is usually referred to Godambe (1960) and McLeish and Small (1988). To 


estimate the SGASC model parameters, we adapt the IFM two-step procedure of Patton (2006) for 
conditional copulas to the MS dynamics. More precisely, the marginals’ conditional distributions 
parameters are estimated in the first step, while the copula GAS parameters are considered in the 
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second step by adapting the Expectation-Maximization algorithm of Dempster et al. (1977). 


The model superior ability to address the issue of tracking the potentially non-linear depen¬ 
dence dynamics usually observed in financial markets is assessed in the empirical part of the work. 
The empirical analysis considers a panel of European regional indexes and focuses on the evaluation 
of the systemic risk contributions of the considered European countries and their evolution during 
the recent GFC of 2007-2008 and the ESDB of 2010-2011. Our analysis confirms that the proposed 
SGASC model is able to explain and predict the systemic risk evolution of the considered countries 
in an effective way. In particular, we find results similar to those recently obtained by pngle et al. 


(2015) using a modified version of the Marginal Expected Shortfall risk measure and Lucas et al. 


(2014b) who consider the European Sovereign debt CDS market. We found that the systemic risk 
contribution of each European country growths during the ESDC of 2010-2011 reaching its highest 
level in mid-2014. After a short decreasing period in the last part of 2014, the magnitude as well as 
the relative importance each country had in terms of systemic risk contribution sharply changed. 
In the very last part of our sample, especially after June 26, 2015 when the Greek government 
unilaterally broke off negotiations with the Eurogroup, we found that the overall systemic risk 
reaches again the mid-2014 levels. 

The remainder of the paper is organised as follows. Section [2] introduces the SGASC model. 
We discuss the marginal models and then we detail the Markov-switching framework to model the 
joint behaviour of the series as well as the dynamic GAS specification. Section [3] concerns the es¬ 
timation methodology and presents the EM algorithm to estimate the GAS parameters. Section [4] 
deals with the systemic risk measurement framework by providing the co-movement risk measures 
for the SGASC model. Section [5] presents data and discusses the main empirical results. Section 
[6] concludes. 


2. The Model 

Let Y t. = (Y\,ti • ■ •, Yd,t ) G be a d-dimensional stochastic vector and let ym-i = (yi,..., yi-i) 
be the past history up to time t — 1 of the weakly stationary stochastic process {Y s ,s > 0}. 
Assume also that exists a first order ergodic Markov chain {S s , s > 0} defined on the discrete space 
Q = {1, 2,..., L} with transition probability matrix Q = {qi t k}, where qi^ = P (St = k \ St -1 = /), 
VI, k G is the probability that state k is visited at time t given that at time t — 1 the chain was 
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in state l, and initial probabilities vector S = (<Ji,..., Sl)' , 81 = P (Si = l ), i.e., the probability of 
being in state l = {1, 2,..., L} at time 1. The SGASC model assumes 


Y t | S t = 0 ~ C (U t ;^ 


(1) 


where U* = ■ ■ ■ t Ud,t) and U^t = Fi (Y),t; &i) denotes the Probability Integral Transfor¬ 

mation (PIT) of Y t f according to its marginal conditional distribution function F) (Y)^; and 
C is the copula distribution conditional to the regime S± = l. Hereafter, we assume 

that generically denotes the parameters of the marginal distribution of Y), for i = 1,2 ,d, 
and the copula distribution may depend on both static G V^p C M m and dynamic k\ G T> k C M n 
state dependent parameters. The SGASC model can be though of as a more general case of the 


static finite mixture of dynamic copulas models of Creal et al. (2013) or an improvement of the 


model of Jondeau and Rockinger (2006). The SGASC specification however differs from those 


classes of models since it introduces a GAS-type dynamic into a finite mixture of copulas, while 
retaining a Markovian structure on the dependence parameters that switches among two or more 


Chollete et al. 

(2009), Pelletier ( 

2006 

) and 

Rodriguez 

(2007 


In this section, we first introduce the general framework to model the univariate marginal con¬ 
ditional distributions and, then, we present the new SGASC specification. It is worth noting that, 
the marginal specification we are going to define should be intended for empirical purposes, and it 
is not strictly related to our new SGASC specification. Indeed, even if the selected marginals’ spec¬ 
ification is highly flexible, and it usually fits well equity returns distribution, alternative marginals 
specifications can be chosen. 

2.1. Marignal models 

We assume each marginal stochastic process s > 0}, i = 1,2,..., d to follows an autoregressive 
process of order 1, AR(1), with Skew Student-t innovations (Fernandez and Steel, 1998) and 


conditional time-varying variance driven by the GJR-GARCH(1,1) specification of Glosten et al 


(1993). Specifically, the AR(l)-GJR-GARCH(l,l)-s5T marginal model can be represented as 


£i,t — 


sST (0, l,Vi,r]i) 


&i,t 

a i,t = + + $2,il(-oo,0] ( £ i,t-l) £ i,t-l + ^3,i a i,t-l^ 
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where (vt, rn) G (0, +oo) x (0, +cx>) represent the degrees-of-freedom and the skewness parameters, 
respectively. The standardised skew Student-t sST (0,1 ,Vi,rji) density for is 


, 0 , 1 j l[0,+oo) ( e i,t) T i tvi {Vi £ i,tt 0 ) 1) 1(—oo,0) ( e i,t) i (5) 

Vi J Vi + ^ 

where £^(-,0,1), denotes the standardised Student-t density with V{ degrees of freedom. To 
preserve the stationarity and to ensure the positiveness of the conditional volatility process a'f f 
for any t = 1,2,... ,T, the following conditions are imposed to the GJR-GARCH(1,1) dynamics 


t Vi (£i,t, 0,1 ,Vi) = 


^ + k 


in equation Wi > 0, j > 0, || < 1, o < d 3 j < 1, with &i ti + d 2 + i?3,i < 1, for i = 


1,2,... ,d, where q = t Vi (e, 0,1, r]i) de = . ^ 2 . The skewing mechanism of Fernandez and Steel 
(1998) provides a flexible way to build a skewed distribution starting from any arbitrary univariate 
unimodal symmetric density function, by including an additional parameter r]i G R + that controls 
for the skewness. The choice of modelling the conditional mean as a first order autoregressive 
process in equation ([2]) is motivated by the need to account for the serial dependence sometimes 


displayed by the financial returns, as described by Embrechts et al. (2003). For easy of readability, 


we group all the marginal model parameters in the vectors Rj = (</>o,j, 4>i,i, $ 2 ,i) Vi j Vi), 


for i = 1,2,..., d. 


2.2. Switching GAS copula model 

Conditionally on the state of the Markov chain St = l, for l = 1,2,..., L, we assume that the 
copula dependence parameter k[ follows the GAS dynamic updating mechanism specified as 

A («t) (6) 

k\ = u l + A z s*_ 1 + B l k l t _ 1 (7) 

= St (k\, V’*) V t (u t ; k[, ^ , (8) 

where u* = (u\^, ■ ■ ■ , u/ G M n and A^, are matrices of parameters of dimension (n x n) 
such that the eigenvalues of B ; are in modulus strictly less than one to preserve the stationarity 
of the GAS dynamic, and is the scaled-score of observation Uj. The parameter ll> 1 controls for 
the state specific level of the GAS dynamic, while the matrices A^ and B ( control for the updating 
step and the persistence of the process, respectively. Moreover, A : M n —> V K in equation ([6]) is an 
absolutely continuous deterministic invertible function that maps M n into the natural parameter 
space V K . When n > 2 it is convenient to specify A as a vector valued function being able to map 
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each component of Kt into the proper space. In general, for the i —th component of A, we consider 
the modified logistic function defined by 


A M (*) = *+ 


, , ( 9 ) 

which maps M into the interval (b, F). The scaled-score function in equation ([8]) can usually be 
obtained as the product of the observation score with respect to the natural parameter K t and a 
proper transformation of the jacobian of the mapping function, in the following way 


S t (k l t ,i/> 1 ) V t (k 1 u iI> 1 ) = z (\ l t ) S t (u t ;K l t,tp l ) V f (u f,K l t ,il > 1 


( 10 ) 


where 


V, 


(u 


dine (u f ;*4, 
dni 


\ = 


dA (k\) 


d k 


and z (•) depends on the choice of scaling mechanism, see 


Creal et al. 


(2013). Hence, Vt (u tj k[, 


denotes the conditional score of the copula probability density function c(-), evaluated at K t and 
Sj is a positive definite, possible parameter-dependent, scaling matrix. Convenient choices 

for the scaling matrix are usually given by 

. 1 —C 


s t ( ) = 




( 11 ) 


where X (kJ, -0*) is the Fisher information matrix that, for well behaved densities, could be written 

d 2 lnc(u t -,K,[,ip l ) 


X[K l t ^ 1 ) =-Et_! 

= E t -i 


Vt X Vt (u 


( 12 ) 


Creal et al. 


(2013) suggest to use the inverse Fisher infor- 


and C is usually set equal to {0, 1}. 

mation matrix, which corresponds to set £ = 1 in equation GD. or its pseudo-inverse square root, 
which corresponds to £ = ^, in order to scale the conditional score for a quantity that accounts for 
its variance. In our empirical tests, we find this latter scaling mechanism much more efficient than 
using an identity scaling matrix (£ = 0). However, sometimes the Fisher information matrix is not 
available in closed form, and we need to resort to simulation or numerical evaluations methods, 
which should be traded-off with approximation degree and code efficiency. In the supplementary 
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material accompanying this paper we report the analytical formulae to compute the score and the 
Fisher information matrix for several bivariate copula specifications: Gaussian, Student-t, Gumbel 
and Clayton. The Frank and Plackett Fisher information matrices are instead evaluated using the 


grid approach proposed by Creal et al. (2013), as detailed in the accompanying supplementary 
material. 


3. Estimation and inference 

As mentioned in the Introduction, model parameters are estimated using a two step procedure that 
consists of a first step where the parameters involved in the marginals are estimated, followed by a 
second step where the dependence parameters are jointly estimated along with the latent Markovian 
states. The resulting Inference Functions for Margins (IFM) two step procedure has been proved to 


be asymptotically consistent for conditional copulas by Patton (2006) to which we refer for further 


details. The IFM procedure can be applied in this context as well as in general dynamic copula 
models since the parameters of the copula distribution in equation ([T]) are separable from those of 
the marginals. In the specific case of the dynamic MS model here considered we also take advan¬ 
tage from the fact that the Markovian dynamic is imposed only on the dependence parameters of 
the copula function and not on the marginals. To estimate the dependence parameters subject to 
the Markovian structure as well as to the GAS dynamic specified in the previous section we adapt 


the Expectation-Maximization algorithm of Dempster et al. (1977). In this section we present the 


EM algorithm, assuming that, the parameters of the marginal distribution have been previously 


consistently estimated by maximum likelihood, see, e.g., Francq and Zakoian (2010). Hereafter 


the pseudo-observations u f are replaced by u t = yF\ , ...,Fd [ydF^djj , where is 

the ML estimate of for all i = 1, 2,..., d. 

The EM algorithm is a powerful and easy programmable tool for ML estimation on data hav¬ 
ing missing structures, such as finite mixtures and Markov-Switching models, and it releases a 
non-decreasing sequence of the log-likelihood function converging to the maximum. For a general 


and up-to-date reference on the EM algorithm see the book of McLachlan and Krishnan (2007). 


In what follows, we present the EM algorithm for estimating the parameters of the SGASC model 
described in Section [2j 

For the purpose of application of the EM algorithm the vector of observations yi-.T, where 
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T is the sample size, is regarded as being incomplete. Following the implementation described in 


McLachlan and Peel 


(2000) the following missing data are consequently introduced zt = zt t 2 ,..., ztjj 

as 


and zzt = (zz t , 1 , 1 , zz t , 1 , 2 , • • •, zz t ^ k ,..., zz ttLiL )' being defined 


zt,l 


ZZ t ,l,k 


1 if St = l, 

0 otherwise 

1 if S t - 1 = l,S t = k , 
0 otherwise. 


Similarly to the latent class approaches, the class membership is unknown and conveniently treated 
as the value taken by a latent Multinomial variable with one trial and L classes, where the temporal 
evolution of class membership is driven by the hidden Markov chain St for t = 1 , 2, ..., T. Aug¬ 
menting the observations yi : r with the latent variables {zt, zzt, t = 1,2, ..., T} allows for replacing 
the log-likelihood function with the complete-data log-likelihood, which becomes 

L L L T 

log£ (5) = Zij log (Si) + EEE zzt,l,k log (qi t k) 

1=1 1=1 k= 1 1 =1 

L T T d 

+ EE ^logc + EE log fi , 

1=1 t= 1 t =1 i =1 


where S = with = (u: 1 , vec(A ; ), vec(B*), i/) 1 ) is a vector containing the parameters of 

the GAS dynamics for the copula dependence parameters n\ and for l = 1, 2,..., L. The EM 
algorithm consists of two major steps, one for expectation (E-step) and one for maximization (M- 


step), see, e.g., McLachlan and Krishnan (2007). On the (m + l)-th iteration the EM algorithm 
proceeds as follows: 


E-step: computes the conditional expectation of the complete data log-likelihood given the observed 
data and the current parameters estimate 

Q(s,H( m )) =E^I^- s(m) )[ log / :( H ) ; yi:T ]. (13) 

M-step: choose Hl m+1 ) by maximizing the preceding expected values with respect to S 

5 (m+1) = arg max Q (a, . (14) 
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The E-step in equation (13) requires the computation of the so-called Q-function, which calculates 


the conditional expectation of the complete-data log-likelihood given the observations and the 
current estimate of the parameter vector Vt = 1,2,... ,T and l = 1,2,... ,L. Exploiting the 
previous factorisation we obtain the following representation of the function Q 

Q (S, 3^) = {log£ (3) ;y 1:T } 

L L L T L T 

OC E h ,1 log (dl) + EEE zz t ,i,k log (qi,k) + ^2 ^2 log 

1=1 1=1 k=l t= 1 1=1 t=l 

where 

Zt,i =F(s t = l\ yi;T) 3^) , zzt,i, k = p (s t -i = l,St = k I yi :T) 3( m )) , 


for l,k = 1,2,..., L, and Vt = 1,2,... ,T, denote the current smoothed probabilities of the states 
evaluated using the well-known Forward-Filtering Backward-Smoothing (FFBS) algorithm de¬ 
tailed in Fruhwirth-Schnatter (2006), and Cappe et al. (2005]). 

The M-step in equation (14) maximizes the function Q (3,3^ m l) with respect to 3 to deter¬ 
mine the next set of parameters E( m+1 ). The updated estimates of the HMM parameters, i.e., the 
vector of initial probabilities S and the transition probability matrix of the hidden Markov chain 
Q are: 

Si=2 zz t,l,k 


fi(m+l) _ 


z l,l 


hm+l) _ 


Ljk= 1 l^t =2 zz t,l,k 

for l,k = 1,2,... ,L, while the parameters 3 can be obtained as the solution of the following 
optimisation problem 

L T 

3 (m+l) _ ar g max EE %logc (u t : S') . 

“ i=i t =i 

Convergence of the algorithm to the ML estimates is guaranteed since the last optimisation step 
delivers a parameter update that increases the log-likelihood function. Standard errors for the 
SGASC dependence and the marginals parameters can be evaluated using the procedure detailed 
in 


Patton (2006), where the Information Matrix of 3 is evaluated numerically after one run of the 


direct Maximum Likelihood estimator for HMM detailed by Zucchini and MacDonald (2009). 
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4. Systemic risk measures 


In this section, we first introduce the two systemic risk measures we consider throughout the paper, 
namely the Conditional Value-at-Risk (CoVaR) and the Conditional Expected Shortfall (CoES) 
and then, we describe how CoVaR and CoES can be calculated assuming the joint returns fol¬ 
low the SGASC model. CoVaR and CoES have been introduced in the systemic risk literature 


by Adrian and Brunnermeier (2011, 2014), and subsequently extended to a parametric dynamic 


framework by Girardi and Ergiin (2013). The CoVaR measures the spillover effects between in¬ 
stitutions by providing information on the Value-at-Risk of an institution or market, conditional 


on another institution’s distress event. In their seminal paper, Adrian and Brunnermeier (2011) 
propose to estimate the CoVaR measure using a system of two quantile equations extending the 


traditional approach of direct Value-at-Risk estimate. Bernardi et al. (2013a, 2015) instead pro¬ 
pose a Bayesian dynamic quantile model where both the VaR and CoVaR equations are function 
of individual and macroeconomic observed risk factors, as in the original CoVaR approach, as well 
as of unobserved components having their own stochastic dynamics. The copula approach to the 


CoVaR has been recently proposed by Reboredo and Ugolini (2015) to evaluate the systemic risk 
in European sovereign debt markets. Here, the CoVaR and CoES are further extended to account 
for both the dynamic evolution of the copula dependence parameters as well as the presence of 
distinct Markovian regimes. The appealing characteristic of the model based CoVaR risk measure¬ 
ment framework is that it inherits the flexibility and easy computability of the dynamic switching 
copula. Hereafter, given the bivariate nature of the CoVaR and the CoES systemic risk measures, 
we consider a bivariate SGASC specification for the stochastic vector (Ym,*) Y),t), j = 1, 2,..., N, 
where “M” denotes the financial System index and “j” denotes the j-th systemically relevant in¬ 
stitution. 

Formally, let (n, 72 ) E [0, l] 2 be predetermined confidence levels, then the CoVaR of the finan¬ 
cial system at time t, denoted by CoVaR^jJ), satisfies the following equation 

IP (*M,t < CoVaR^jJ 2 | Y j>t < VaR} 2 ) = n, (15) 

where VaR)" 2 denotes the marginal Value-at-Risk (VaR) of institution j such that P (Yjj < VaR)) 2 ) = 


r 2 . Roughly speaking, the Conditional Value-at-Risk of the financial system, is the quantile of the 


distribution of conditional on an extreme event affecting institution j’s returns Yj,t- As in Gi- 
























rardi and Ergiin (2013), we define such an extreme event as Yjt being below its VaR at confidence 
level 72- 


Remark 4.1. The definition of Conditional Value-at-Risk in equation (15) is substantially dif¬ 


ferent from that originally presented in Adrian and Brunnermeier (2011) and coincides with that 


proposed in Girardi and Ergiin (2013). As discussed in Girardi and Ergiin (2013) and Mainik and 


Schaanning (2014) this definition essentially preserves the stochastic ordering introduced by the 


joint distribution. 


Given the dynamic context introduced in the previous Sections, the random variable we refer to for 
the calculation of the forward looking systemic risk measure CoVaRjljJ 2 +1 is (Yjyi,t+i , Yjj+i I ynt, St) 
The CoVaRj^j r) +1 and VaRj 2 +1 characterise the conditional and marginal quantiles of the predic¬ 
tive distribution of (Vm^+Ii ^j,t+ 1 I ynt, St). The next Proposition provides the natural link between 
the dynamic model and the approach to systemic risk assessment by characterising the predictive 
distribution. 


Proposition 4.1. LetYt = (4m,«+1) Yj,t+\) follows a bivariate SGASC process, then the one-step 
ahead predictive cumulative distribution function ofYt at time i + 1 , given information up to time 
t is a mixture of component specific predictive cumulative distributions 

L 

H (Y m | y 1:t , S t ) = Y *th\t C (Un-i I S t+ 1 = h yi a ), (16) 

i=i 

with mixing weights 

L 

^t+i\t = Y q m,i F (S t = m I yi:t), l = 1,2 ,... ,L, (17) 

771=1 


where U t+ i = (F M (V M ,t+i, ^m) , Fj {Y jtt+ i,&j)), j = 1,..., N and q m j, is the (m, l)-th entry of the 


Markovian transition matrix Q, see, e.g., Zucchini and MacDonald (2009). 


It is worth adding that equation (16) follows immediately from the fact that the cumulative pre¬ 
dictive distribution of Y t+ i given the past history of the process is a finite mixture of copulas and 


that mixtures of copulas are copulas themselves, see, e.g., Durante and Sempi (2015). We now 
provide expressions for the one-step-ahead “predictive” systemic risk measures. 

In principle, the calculation of CoVaR requires the prior evaluation of institution’s j marginal 
VaR that can be performed by inverting the marginal cdf of Yj, i.e., VaRj 2 t+1 = Ffi 1 ( 72 ; Rj). From 
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a computational point of view, the copula approach is more tractable because it does not require 
the evaluation of the marginal VaR which coincides with the PIT at the chosen quantile confidence 


level 72- Conditional on VaRj 2 +1 , the CoVaR 




is calculated as the value of 


Vt+i 


such that 


(V M ,m < yl +1 ,Y jtt+ 1 < VaRJ 2 t+1 ) = nr 2 , 


(18) 


or in terms of the bivariate SGASC predictive distribution we consider here 
L 


E 

l=i 


J l ) C 

7T t+l\t L ' 




(em (y*+i,^M) 

L 

= Y I\t c ( Fm (y*+ 1 ) ^m) , r 2 ; K l t, = TrT- 2 . 


(19) 


/=! 


Adrian and Brunnermeier (2014) also propose to extend the expected shortfall (ES) risk measure 
to a systemic framework by evaluating the marginal ES of each institution at the CoVaR level. 
The CoES is defined as the expected shortfall of Ri.t+i below its CoVaR^|J 2 +1 level, conditional 
to Y],t+\ being below its VaRj 2 +1 level, i.e., 


CoES ™ +1 


= ES (y M> t+i < CoVaR^jg 2 +1 I Y ht+1 ^ 


< VaRj 2 +1 


= E (y M ,t+i | E M , m < CoVaR^jJ 2 +1 ,y i)t+ i < VaRj 2 t+1 


( 20 ) 

( 21 ) 


Our definition of CoES coincides with the one given by Bernardi et al. (2013b), Bernardi and 


Petrella (2015) and Mainik and Schaanning (2014) and substantially differs from the one of Adrian 


and Brunnermeier 


(2011). Within the SGASC framework, the forward looking CoES^jJ 2 +1 can be 


evaluated by numerical integration of the CoVaRllT 2 „ as follows: 


Oo E SSiS +1 = ^/" C »VaRSiJ, t+1 


d'y. 


( 22 ) 


As discussed by Bernardi et al. (2013b) and Mainik and Schaanning (2014), the CoES risk measures 


inherits the same properties of the Expected Shortfall (ES) such as the sub-additivity with respect 
to linear combinations, see, e.g., Artzner et al. (1999 ). As a direct consequence of the sub-additivity 
property, the CoES can be effectively used in order to measure the total systemic risk contribution 


of different assets to the overall financial system. In a different context, Engle et al. (2015) suggest 


to employ a linear combination of individual Long Run Marginal Expected Shortfall (LRMES) in 
order to obtain an aggregate measure of the total market systemic risk for the European region. 
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Using same arguments, Brownlees and Engle (2015) also rely on the sub-additivity property of the 
LRMES to get an aggregate version of the Systemic Risk (SRISK) indicator. From our definition 


of CoES in equations (20) (21), it is easy to see that the LRMES can be obtained as a special 
case of the CoES by simply letting CoVaRj))jJ 2 +1 —> oo, i.e., by imposing T\ = 1. To aggregate 
the individual systemic risk levels to get an overall indicator of total systemic risk for the whole 
economy, we can define the total market forward looking CoES as 

N 

„„ .r^T?c 


C°ES™ ,=£w,CoES™ 


(23) 


3 =1 


where the scalars Wj, j = 1, 2,..., N, denote the weights associated to each one of the N systemi- 
cally relevant financial institutions belonging to the market. The definition of CoES^’J^ is useful 
to estimate the total loss the overall market is going to face as a consequence of a crisis affecting a 
market participant, which is transmitted to the market when we observe a realisation of 1 m below 
the CoVaR level. 


As discussed in Adrian and Brunnermeier (2011, 2014), Girardi and Ergiin (2013) and Mainik 


and Schaanning (2014), it is also useful to consider the difference between the CoVaR and the CoES 
from their “median” value. Here, we consider the ACoVaR^jJ 2 t+1 and the ACoES^jJ 2 +1 quantities 
defined as 


ACoVaR£j5 +1 = 100 x 


CoVaR: 


n,T2 


- CoVaR 


ri,b J 


CoVaR: 


Ti,b J 

M\j,t+1 


and 


ACoES M|S + i = 100 X 


CoES 


MRt+l 


— CoES 


T1,V> 


c°Es^; m 


(24) 


(25) 


where V represents the benchmark state that we define as P (lyt+i < VaR^f +1 ) = 0.5, i.e., the 
CoVaR of the system when the country specific indexes are below their median value. The ACoVaR 
and the ACoES measure the percentage increase of the systemic risk conditional on a pre-specified 
distress event. It follows that the ACoVaR and the ACoES can be effectively used to measure 
how the CoVaR and the CoES change when a particular institution becomes financially distressed. 
In other words, the ACoVaR and the ACoES estimate the dynamic evolution of the specific 
institution j’s contribution to the overall systemic risk. Furthermore, these two quantities can be 


employed for policy rules and for risk management. Adrian and Brunnermeier (2011, 2014) found 


strong evidence for the existence of a relation between the ACoVaR and several macroeconomic 
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indicators. In the empirical part of the paper, we will employ the ACoVaR and the ACoES risk 
measures to investigate the systemic risk contributions of specific countries to the overall risk of 
the European economic system. 


5. Empirical study 


In what follows, we apply the econometric framework and the methodology described in the pre¬ 
vious sections to examine the evolution of the systemic risk in Europe over the past decade. As 


discussed in Bernardi et al. (2013bl), Bernardi and Petrella (2015) and Billio et al. (2012) among 


others, the high level of interdependence and interconnection among financial institutions has been 
recognised as the main ingredient that facilitates adverse shocks affecting individual institutions or 
countries to spread over the overall financial system. Systemic events are particularly relevant since 
they involve extreme losses for all the market participants threatening the stability of the entire 
economic and financial system. The evaluation of systemic risk in the European financial system 


has been recently considered by Reboredo and Ugolini (2015) in order to assess the impact of the 


recent ESDC of 2010-2011. Billio et al. (2012) analyse stock price data on hedge funds, banks, 


brokers and insurers for both the US and the Euro market. More recently, Engle et al. (2015) and 


Lucas et al. (2014b) analyse the evolution of the European systemic risk using a new systemic risk 
indicator (LRMES) and the multivariate Generalised-Hyperbolic GAS model, respectively. 


5.1. Data 

To investigate the systemic risk contribution of several European countries to the overall European 
system we consider eleven equally weighed portfolios composed by the fifteen most capitalised 
companies domiciled in each country. By construction, the considered country specific indexes 
are representative of the equity market of the country they belong to. The selected countries 
are Austria (AU), Belgium (BEL), Denmark (DEN), France (FRA), Germany (GER), Hungary 
(HUN), Italy (IT), Netherland (NET), Spain (SPA), Sweden (SWE) and United Kingdom (UK). 
To evaluate the systemic risk we also need to identify a proxy for the total Eurozone equity market 
(MKT) which consists of an equally weighed portfolio across all the companies belonging to our 
dataset. Details about the composition of the regional indexes are reported in Table [A~2 

We analyse equity indexes’ log-returns from July 8, 1999 to October 16, 2015, covering the 
recent Global Financial Crisis of 2007/2008 as well as the European sovereign debt crisis of 2010. 
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For each index, the last H = 2000 observations, covering the period from November 12, 2007 to 
the end of the sample, are considered to perform the out-of-sample forecasting exercise, while the 
first part is used to estimate models’ parameters and to assess models’ performances with respect 
to nested alternatives. Descriptive statistics for all the considered series are reported in Table 
A.3. In line with most important stylised facts, frequently detected in financial time series, the 
returns appear to be negatively skewed and leptokurtic, indicating that their empirical distributions 
strongly deviate from the Gaussian one. The departure from normality is confirmed by the Jarque- 
Bera (JB) statistic which always rejects the null hypothesis at the 1% level of significance. The 
presence of large volatility clusters followed by periods of low volatility is also well documented 
by the data. These facts are coherent with the presence of different regimes of “bull” and “bear” 
market conditions, which are usually associated with “high” and “low” dependence, as discussed 
in 


Pelletier (2006) among others. 


5.2. Marginal specifications 

In order to filter out all the stylised facts affecting the univariate marginal conditional distributions 


we estimate the AR(l)-GJR-GARCH(l,l)-s<ST" detailed in Section 2.1 Estimated coefficients of 
each marginal model are reported in Table [A~4 Our results seems to be in line with those usually 
found in the financial econometric literature, such as the strong persistence and the positive reaction 
of the conditional volatility to past negative innovations. The skewness parameters fg of the Skew 
Student-t distribution, are always significantly smaller then one, justifying our choice of skew 
innovations. The estimated degree of freedom parameters Vi strongly confirms the excess of kurtosis 
and the departure from normality of each of the considered series. To check the goodness of fit 
of the estimated marginal distributions, we test if the PITs implied by the estimated conditional 
densities are independently and identically distributed uniformly in the unit interval (0,1). To this 


end, we implement the same testing procedure employed by Vlaar and Palm (1993), Jondeau and 


Rockinger (2006) and Diebold et al. (1998). Specifically, the iid uniform test is made of two main 


parts. The first part checks the independence assumption by testing if all the conditional moments 
of the data up to the fourth one have been captured by the model, while the second part aims 
to verify whether the Skew Student-t assumption is reliable by applying a Uniform (0,1) test on 


the PITs. The first test consists to examine the serial correlation of the quantities (ui,t ~ Uij k for 
k = 1, 2,..., 4, where Ui = T^ 1 Ylt =i by performing a simple linear regression of (u^t — Ui) k 
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on their 20 lags. The null hypothesis of absence of serial correlation is tested using a Lagrange 
multiplier-type test defined by the statistics (T — 20) R| where R|, with k = 1,2, ...,4 is the 
coefficient of determination of the regressions. The corresponding test statistics DGT-AR(/c) for 


k = 1,2, 3,4, are reported in the first three rows of Table A.5 We note that, for almost all the 
considered series, the test is in favour of the null of absence of serial dependence for the first four 
conditional moments of the estimated PITs. Concerning the uniform distribution test, following 


Jondeau and Rockinger (2006), the statistic tests have been evaluated by splitting the empirical 


distributions in G = 20 bins. For more information about the test implementation, see, e.g., 
Diebold et al.| (T998|) and Jondeau and Rockinger (2006). The estimated statistics for the uniform 


test, DGT-H(20), are reported in the last row of Table A.5 Except for DEN, HUN, SWE and 
UK, for the remaining country indexes the test indicates that the estimated PITs are uniformly 
distributed over the interval (0,1) at a confidence level lower then 1%. Concerning HUN and SWE, 
the rejection of the uniform assumption for the PITs is principally related to the large number of 
zeros those series exhibits, especially in the very first part of our sample. This empirical evidence 


is also displayed in Figure A.l where the empirical distribution of the PIT series along with the 
5% approximated confidence levels are reported. However, since this finding is quite common and 
only affects the centre of the marginal distribution, we decide to not pre-filter the series and to 
continue with our empirical investigation. 

5.3. Dynamic copula specifications 

The SGASC model nests several alternative copula specifications: it reduces to the GAS copula 


(GAScop) model of Creal et al. (2013) and Harvey (2013) when L = 1, while for a 1 = fi 1 = 0, 
V7 = 1,2, ...,L, it reduces to the Markov Switching copula model of Jon deau and Rock inger 


(2006) and 


Chollete et al. 


(2009) (MScop), and to the static copula (STATcop), for a 1 = (3 l = 0 


and L = 1. Concerning the choice of the copula distribution, in an unreported analysis, we find 
that the Student-t copula seems to be the best choice in order to describe the dependence patterns 
across financial indexes. This evidence is quite common in the financial econometric literature, see, 


e.g., Jondeau and Rockinger (2006), Rodriguez (2007) and Demarta and McNeil (2005). To select 


the best model, we compare the GAScop, the STATcop as well as the MScop and the SGASC 
models with different number of regimes, namely L = 2,3,4. For each pairs of country and market 


indexes, Table A.6 reports the Akaike Information Criteria (AIC) for all the Student-t copula 
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specifications we consider. The information criteria select the SGASC-L2 specification for 6 coun¬ 
tries, while the simpler GAScop model is preferred in the remaining five cases. This evidence 
suggests that countries like France, Italy, Netherlands do not display abrupt changes in their de¬ 
pendence structure with the overall European equity index while Markovian patterns are detected 
for countries such as, for example, Belgium, Denmark and Germany. Table |A.6| reveals also that 
the standard static copula specification STATcop is clearly suboptimal when compared with more 
flexible models. 

Before turning to the out-of-sample forecast exercise, two additional points should be ad¬ 
dressed. The first concerns the goodness of fit analysis of the copula specification, while the latter 
refers to the adequacy of the Markov switching dynamics. As regards the goodness-of-fit analysis, 
we aim to tests the null hypothesis of correct copula specification 


H 0 : C (U ht , U 2 , t ) = C 0 (U ltt , U 2 ,t ), Vt = 1,2,. 


T 

5 1 5 


(26) 


with respect to the true unknown copula distribution Cq ( U\t , U 2 t)- Table A.7 reports the p-values 


of the Anderson-Darling test statistics performed using 1000 parametric bootstrap resamples as 


indicated in Manner and Reznikova (2012) and Patton (2012), to which we refer for further details 


Not surprisingly, except for Belgium, the goodness-of-fit test results are in line with those obtained 
by the model selection criterium AIC. This evidence supports our choice of performing the systemic 
risk analysis using those models selected by the AIC. The latter analysis investigates which SGASC 
parameter drives the Markov-Switching behaviour of the dependence structure. Specifically, we test 
the hypothesis of regime independence for the parameters {oj, a, f3, is} using a standard likelihood 


ratio test. LR p-values are reported in Table A.8 Once again the LR results strengthen the 
evidence provided by both the AIC and the goodness-of-fit test. At both extremes, there is no 
evidence of state dependence for Hungary, while United Kingdom displays well identified regime 
dynamics. In between those cases, we find indexes with low or moderate Markovian structure. 
For example, Belgium and Germany display significant switches in both the unconditional mean 
and persistence of the dependence parameter, while Denmark is characterised by different tail 
regimes since only the degrees-of-freedom parameter changes according to the regime. Finally, 
Table A.9| reports the in-sample estimated coefficients for the selected models. Concerning the 
SGASC specifications, we observe that almost all coefficients are strongly significant with high 
persistence in each state. The positive and significantly different from zero impact of the scaled 
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scores to the copula parameters, suggests that the GAS dynamic effectively moves the dependence 
parameters to the proper direction. Those findings hold also for the simpler GAScop specification. 


5-4- Systemic risk contributions 

In this section, we apply the estimated models and the systemic risk measures introduced in 
Section [4] to asses the systemic risk contribution of each country index to the overall European 
equity market. To this end, one step ahead rolling forecasts are performed using H = 2000 out- 
of-sample observations covering the period beginning on November 12, 2007 to the end of the 
sample. Parameter estimates are updated every 25 observations for a total of 80 refits. The 
estimated coefficients are quite stable during the forecast period suggesting that the underline 
structure of the economy is well described by the selected models. The estimated coefficients over 
the out-of-sample period are not reported to save space and are available upon request to the 
second author. The selected models are employed to predict the CoVaR and the CoES of the 
overall European system M, measured by the Market index (MKT), conditional to distress events 
affecting each European region, measured by the corresponding country indexes, i.e., CoVaR^Jj 2 
and CoES^jj 2 , for j € {AU, BEL, DEN, FR, GER, HUN, IT, NET, SPA, SWE, UK}. The confidence 
levels (ti,T 2 ) are fixed at 5% which means that we use the situation where the country specific 


index is below its 5% marginal Value-at—Risk level as conditional distress event. Figure A.2 
reports, for the forecast period, the predicted CoVaR^jj 2 and CoES^jj 2 . Vertical dashed lines 
represent the major financial downturns experienced by the European economic system during the 
period 2007-2014. A timeline of the major European financial crisis is provided in Table |A.1| 
In the bottom panel of each subfigure, for the SGASC specifications, we report the estimated 


out-of-sample smoothed probabilities, ¥ (St = 1 | Yi : t+h)- Figure A.2 gives insights about the 
dynamic evolution of the systemic risk contributions during the different economic and financial 
phases the European system experienced since the end of 2007. As documented by Figure [A~2j the 
CoVaR and the CoES systemic risk measures suddenly adapt to the relevant changes affecting the 
underlying financial system’s wealth, such as those experienced during the turbulent phases before 
the bails out of Portugal and Greece in 2011. Concerning instead the interpretation of the smoothed 
probabilities of the hidden Markov chain in Figure pV2| it is worth stressing that the dynamic non¬ 
linear evolution of the copula dependence parameter as well as that of the degrees-of-freedom 
prevents a clear and straightforward identification of the latent states with given characteristics 
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of the conditional distribution. To further explain the nature of the identification problem, let 
consider the possibility that for some assets we could observe, for example, that high levels of linear 
dependence are associated with high values of the degrees-of-freedom parameter. In those cases, 


the interpretation of low and high probability levels in the bottom panels of Figure A.2 becomes 
cumbersome. Furthermore, another relevant point to be addressed is that the interpretation of 
the states is not homogeneous across countries. As an example, Denmark in Figure [A . 2c| suddenly 
experienced an abrupt decrease of the smoothed probability from a state of low correlation/high 
tail dependence to a state of high correlation/moderate tail dependence. This is not the case, for 
example, of Germany where, during the same period, the system switched from a state of high 
correlation/high tail dependence to a state of low correlation and low tail dependence. To gain 


further insights about the identification issues, we report in Table A. 10 the long run dependence 
measures as well as two measures of persistence in terms of expected duration, conditional to the 
Markovian state. Concerning the conditional expected durations, we report the half life (HL) index 
as a measure of speed reversion of the dependence parameter dynamics and the implied duration 
of each regime (D) of the Markov chain. Moreover, we also consider the long run correlation pi and 
tail dependence zui, for l = 1,2. Except for Austria, which is characterised by moderate long run 
tail dependence in each regime, all the remaining SGASC specifications clearly identify two regimes 
of high and low long run tail dependence and correlation. European countries seams to be also 
quite heterogeneous with respect to the persistence of the conditional dependence parameter and 
the expected duration of the Markov chian. Indeed, Spain and Germany display high persistence 
in the conditional dependence dynamic but low persistence in the Markov chain. On the contrary, 
Austria reports low persistence in the conditional correlation but very high persistence in the 
Markov chain. These findings suggest that different European countries react differently to the 
new information coming from the market. Moreover, we note that SGASC models usually display 
higher persistence in terms of HL of the conditional process compared to the simpler GAScop 
specifications. Indeed, the additional flexibility gained from the switching behaviour of SGASC 
models, permits to disentangle the information that impacts only on the memory of the dependence 
process and that impacting only on some characteristic of the dependence structure such as the 
levels of correlation and tail dependence. 
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5.5. Backtesting the Co VaR systemic risk measure 

Now we move to the evaluation of the CoVaR forecasts. Since the CoVaR risk measure is essentially 
a modified version of the Value-at-Risk, we can employ the usual VaR backtesting procedures such 


as the unconditional (UC) and conditional coverage (CC) tests of Kupiec (1995) and Christoffersen 


(1998]), the Actual over Expected (AE) ratio and the mean and maximum absolute deviation 


(ADrnean, ADMax) considered by McAleer and da Veiga (2008). The only difference between our 
context and the usual VaR backtesting procedure, concerns the evaluation of the “hitting sequence” 


of returns exceeding the CoVaR levels. In line with our definition of CoVaR in equation (15), we 


follow the approach of Girardi and Ergiin (2013) who consider the series of returns jointly exceeding 


both the CoVaR and VaR levels as the proper hitting sequences. In Table [X. 11 1 we report, for each 
country, the AE, ADrnean, ADmax, as well as the p-values of the UC and CC tests. We observe 
that, except for BEL, the null hypothesis of correct coverage of the conditional lower tail of the 
joint distribution at the 5% confidence level is never rejected. ADrnean suggests that, on average, 
CoVaR violations are of the order of magnitude of about 0.3%, while the ADmax suggests that, 
during period of distress, CoVaR absolute deviations could exceed the 3.5%. 

5.6. Ranking European countries 

The ACoVaR and ACoES risk measures are particularly useful in order to rank the countries in 
terms of their systemic risk contributions. Note that, the resulting ranking does not tell anything 
about the intrinsic riskiness profile of the particular country, but instead, it is informative about 
the role each country plays in the overall European system. This point is particularly relevant 
when systemic risk analyses are carried on. Indeed, the countries on the top of the rank are going 
to be those who play a central role into the overall European equity market. The daily evolution 


of the ACoVaR and ACoES, calculated as in equation (24) and (25), are plotted in Figure A.3 


along with the superimposed smoothed estimates. The visual inspection of the figure reveals three 
relevant phenomena. First, as expected, both risk measures provide the same underlying systemic 
risk signal, with the ACoES being always below the ACoES with almost the same dynamic pattern. 
Second, ACoVaR and ACoES seem to display a quite heterogeneous behaviour across countries. 
Third, the systemic risk importance of the considered countries increases during the second part of 
the sample. The only two exceptions are Austria, which has constant systemic risk contributions 
over the whole sample, and Hungary, which instead is less systemically important over the period 
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2011-2015. Figure [A~5l reports the dynamic evolution of the out-of-sample monthly averages of the 
ACoVaR and the ACoES risk measures, and provides insights on the relative importance of each 
country to the overall European equity market. Monthly ACoVaR and ACoES are calculated by 
averaging the daily levels plotted in Figure |A.3[ Concerning HUN, the ACoVaR and the ACoES 
risk measures are substantially lower than the estimated levels of others countries and are not 


reported in Figure A.5 According to the ranking, over the last decade, France is the country with 
the highest systemic risk contribution, followed by Spain, Germany and Netherlands. For those 
countries the ranking remains unchanged over the whole out-of-sample period, while for other 
countries such as, Italy and UK we observe a huge increase and decrease of the ACoVaR just 
after the onset of the European crisis, respectively. During the last part of the sample we observe 
an increase of the systemic risk even for AU and DEN, which were at the bottom of the ranking 
prior to the European crisis. Our results are strictly in line with the findings recently reported by 


Engle et al. (2015). To summarise, we found that the systemic risk contribution of each European 


country grows during the ESDC of 2010-2011 reaching its highest level in mid-2014. After a short 
decreasing period in the last part of 2014, the magnitude as well as the relative importance each 
country has in terms of systemic risk contribution, sharply changed. In the very last part of our 
sample, just after the decision of the Greek government to unilaterally broke off negotiations with 
the Eurogroup on June 26, 2015, the ACoVaR and ACoES measures suddenly rise up to the highest 
levels experienced since the mid-2014, for all the European countries. It is worth mentioning that 
the main difference between the ACoVaR and the ACoES concerns the end of the out-of-sample 
period. Indeed, during this subperiod the systemic signal provided by the ACoES risk measure is 
clearly better then that of the ACoVaR which displays the tendency to shrink toward the same 
level. 


5.7. Dependence measures forecast and contagion effects 

The SGASC model can be effectively used to forecast several measures of association between the 
random variables such as the linear correlation, the concordance measures Spearman’s rho and 
Kendall’s tau and the coefficients of tail dependence. The Markovian nature of the SGASC model, 
and its flexibility to accommodate the marginals’ specification, entails that the evaluation of such 
dependence measures can become cumbersome. The only exception concerns the evaluation of 
upper and lower tail dependence. In fact, since these measures only depend on the copula specifi- 
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cation, the tail dependence coefficients of the predictive distribution for SGASC models are going 
to be equal to the convex linear combination of the components specific ones, conditional to the 
Markovian states. The linear correlation coefficient, the Spearman’s rho and the Kendall’s tau, 
instead, are not analytically available, and the only achievable solution is to rely on simulation 


methods. Following the approach of Chollete et al. (2009), we evaluate the linear correlation, the 


Spearman’s rho and the Kendall’s tau coefficients empirically based on 10000 simulated draws from 


the joint predictive distribution defined in equation (16) and we repeated the procedure for the 


whole out-of-sample period. Figure |A.4 reports the predicted dependence measures, along with 


the corresponding filtered marginal volatilities. Figure A.4 provides clear evidence of changing 
dependence structures during periods of financial turmoil, revealing also that the SGASC model 
adequately predicts the upward shifts in the dependence measures at economically relevant dates 
represented by vertical dashed lines. In this respect, particularly evident is the onset of the Eu¬ 
ropean sovereign debt crisis of 2010 which coincides with a suddenly increase of the dependence 
measures for all the investigated countries. Furthermore, by comparing Figure |A.4 and A.5 it 
emerges that the ordering induced by the lower tail dependence measures almost corresponds to 
that introduced by the systemic risk contributions. The only exception concerns UK which reports 
high tail dependence and moderate levels of both ACoVaR and ACoES. A possible explanation 
for this result should be ascribed to the fact that tail dependence at penultimate levels may be 


significantly stronger than in the limit, see Manner and Segers (2011). Consistently with previous 


results, Hungary reports a close to zero tail dependence index over the whole out-of-sample period. 


6. Conclusion 

In this paper, we propose a new MS dynamic copula model (SGASC) that explicitly accounts for 
the presence of different Markovian regimes as well as a smooth within-regime dynamic evolution of 
the dependence parameters. Specifically, exploiting the recent advantages for score driven processes 


of Creal et al. (2013) and Harvey (2013), we allow the state dependent copula parameters to be 


updated using the scaled score of the conditional copula distribution. This choice allows us to 
achieve a greater level of flexibility in the context of dynamic copula models, and it also introduces 
a stochastic behaviour for the class of GAS models in a natural and effective way. The SGASC 
model nests several alternative copula specifications with increasing levels of complexity such as 
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the static copula, the Markov switching copula of Jondeau and Rockinger (2006) and Chollete et al. 


(2009) as well as the GAS copula specification of Creal et al. (2013) and Harvey ( 2013| ). Model 
parameters are consistently estimated by employing the Inference Function for Margins (IFM) two 


step estimator, where the second step is performed by tailoring the EM algorithm of Dempster et al. 


(1977) to this class of models. The proposed estimation methodology is consistently equivalent to 
a single step estimation as long as the marginal conditional distributions do not depend on the 
latent Markovian states. 

Although the SGASC model can be effectively used to understand the dynamic evolution of 
the non linear dependence among financial assets, here we focus our attention on the systemic risk 
measurement. Indeed, another relevant contribution of the paper is to introduce and estimate the 
Conditional Value-at-Risk (CoVaR) and the Conditional Expected Shortfall (CoES) risk measures 
of Adrian and Brunnermeier| (2011, 2014) and Girardi and Ergiin (2013), for the class of SGASC 
models. 

In the empirical part of the paper, we provide a comprehensive study to analyse the evolution of 
the systemic risk in Europe over the past decade. The considered period includes the recent Global 
Financial Crisis 2007-2008 and the European Sovereign Debt Crisis of 2010-2011. As regards our 
results, we found empirical evidence in support of the SGASC specification for more than half of 
the selected European countries, while the remaining countries are found to be well represented 


by the simpler GAS copula specification of Creal et al. (2013) and Harvey (2013). The country 
individual systemic risk contribution as well as the overall systemic risk level in Europe is evaluated 
by means of the CoVaR and CoES risk measures and their A counterparts. Our empirical results 
confirm that, the proposed SGASC model is able to explain and predict the systemic risk evolution 


in an effective way. According to Engle et al. (2015), France is found to be the most systemically 
important country followed by Spain, Netherlands and Sweden. We find that the systemic risk 
contributions grow during the ESDC of 2010-2011 reaching its highest level in mid-2014. Then, 
after a short decreasing period in the last part of 2014, the magnitude as well as the relative 
importance of each country’s systemic risk contribution, sharply changed. In the very last part of 
our sample, just after the decision of the Greek government to unilaterally broke off negotiations 
with the Eurogroup on June 26, 2015, the ACoVaR and ACoES measures suddenly rise up to the 
highest levels experienced since the mid-2014, for all the European countries. Moreover, we find 
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evidence that, during this final period, the ACoES risk measure provides a more clear systemic 
ranking. To conclude our empirical analysis, we also investigate the evolution of several dependence 
measures. Our main finding reveals that the ordering of systemic importance induced by the 
ACoVaR and ACoES risk measures is almost preserved by the lower tail dependence coefficient. 
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Appendix A. Tables and figures 


Date Event 

January 21, 2008 the global stock markets suffer their largest fall since September 2001. 

March 16, 2008 the Bear Stearns acquisition by JP Morgan Chase. 

September 15, 2008 the Lehman’s failure. 

March 9, 2009 the peak of the onset of the recent GFC. 

December 8, 2009 the downgrading of Greece’s credit rating from A- to BBB+ by Fitch ratings agency. 

May 18, 2010 the Greece achievement of 18bn USD bailout from EFSF, IMF and bilateral loans. 

November 29, 2010 the Ireland achievement of 113bn USD bailout from EU, IMF and EFSF. 

May 05, 2011 the Portugal bailout from ECM. 

August 05, 2011 the S&P downgrading of US sovereign debt. 

March 16, 2013 the Cyprus achievement of 13bn USD bailout from ECM. 

April 07, 2013 the conference of the Portuguese Prime Minister regarding the high court’s block of austerity plans. 

April 30, 2013 the approval of the Cyprus bailout by the Euro Parliament. 

September 17, 2013 the drop of car sales to the lowest recorded level in the Euro area. 

June 03, 2014 the drop of Eurozone inflation, and the consequent increasing pressure on the Central Bank. 

September 09, 2014 the mediterranean countries prepare for further unrest. 

November 28, 2014 the Italian unemployment rate reaches the record high since the 1977. 


June 26, 2015 the Greek government unilaterally broke off negotiations with the Eurogroup. 


Table A.l: Financial crisis timeline. 
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Name 

Min 

Max 

Mean 

Std. Dev. 

Skewness 

Kurtosis 

1% Str. Lev. 

JB 

Kendall’s r 

p(i) 

In-sample, from 

08/07/1997 to 09/11/2007 







Austria 

-6.19 

3.47 

0.05 

0.69 

-0.91 

10.37 

-1.88 

5101.1 

0.35 

0.02 

Belgium 

-3.98 

6.33 

0.02 

0.88 

0.1 

7.17 

-2.53 

1542.75 

0.51 

0.11 

Denmark 

-6.71 

4.86 

0.06 

1.08 

-0.3 

5.32 

-2.8 

508.47 

0.5 

0.02 

France 

-7.63 

7.12 

0.01 

1.27 

-0.17 

6.14 

-3.66 

882.71 

0.67 

0.05 

Germany 

-5.34 

6.19 

0.02 

1.22 

-0.11 

5.47 

-3.36 

542.83 

0.66 

-0.01 

Hungary 

-4.77 

4.08 

0.05 

1 

-0.14 

4.27 

-2.62 

149.62 

0.34 

0.03 

Italy 

-7.75 

7.48 

0.01 

1.18 

-0.2 

7.09 

-3.53 

1499.85 

0.65 

0.03 

Netherlands 

-5.31 

5.93 

0.01 

1.24 

-0.07 

4.93 

-3.36 

331.75 

0.67 

0.02 

Spain 

-7.2 

6.09 

0.05 

1.03 

-0.35 

5.92 

-2.79 

797.23 

0.59 

-0.01 

Sweden 

-7.01 

6.21 

0.04 

1.24 

-0.15 

5.86 

-3.35 

730.84 

0.64 

0.01 

United Kingdom 

-6.57 

5.85 

0.02 

1.18 

-0.22 

5.95 

-3.58 

786.92 

0.62 

0 

Market 

-5.13 

4.61 

0.03 

0.86 

-0.38 

6.12 

-2.57 

910.75 

- 

0.05 

Name 

Min 

Max 

Mean 

Std. Dev. 

Skewness 

Kurtosis 

1% Str. Lev. 

JB 

Kendall’s r 

p( i) 

Out-of-sample, from 09/11/2007 

to 16/10/2015 







Austria 

-7.15 

6.77 

-0.01 

1.42 

-0.37 

6.03 

-4.33 

812.77 

0.65 

0.12 

Belgium 

-5.4 

6.14 

0.01 

1.2 

-0.11 

5.69 

-3.55 

610.77 

0.73 

0.09 

Denmark 

-9.49 

8.44 

0.02 

1.47 

-0.22 

7.56 

-4.41 

1756.1 

0.64 

0.07 

France 

-9.32 

10.1 

0 

1.63 

-0.17 

6.88 

-4.99 

1266.06 

0.79 

0.05 

Germany 

-7.62 

9.96 

0 

1.42 

-0.03 

8.2 

-4.48 

2259.45 

0.73 

0.02 

Hungary 

-8.83 

6.6 

-0.03 

1.06 

-0.63 

11.3 

-3.22 

5886.51 

0.39 

0.06 

Italy 

-7.48 

8.81 

-0.03 

1.83 

-0.17 

4.5 

-5.09 

198.35 

0.71 

0.01 

Netherlands 

-7.04 

7.3 

0.02 

1.37 

-0.29 

7.01 

-4.09 

1372.06 

0.74 

0.02 

Spain 

-7.82 

11.03 

-0.02 

1.71 

-0.07 

5.91 

-4.74 

709.44 

0.71 

0.06 

Sweden 

-7.76 

8.9 

0.02 

1.46 

-0.03 

7.49 

-4.48 

1687.4 

0.7 

-0.01 

United Kingdom 

-8.83 

6.93 

0.01 

1.36 

-0.22 

6.73 

-3.84 

1178.19 

0.69 

0 

Market 

-6.9 

7.6 

0 

1.29 

-0.24 

6.93 

-3.89 

1311.09 

- 

0.06 


TABLE A.3: Summary statistics of the panel of country specific indexes along with the total Market, for the period 
beginning on July 8, 1999 and ending on October 16, 2015. The seventh column, denoted by “1% Str. Lev.” is 
the 1% empirical quantile of the returns distribution, while the eight column, denoted by “JB” is the value of the 
Jarque-Bera test-statistics. The last two columns report the estimated Kendall r with respect to the total Market 
and the first order empirical autocorrelation of returns, respectively. 
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Austria 

Belgium 

Denmark 

France 

Germany 

Hungary 

Italy 

Netherlands 

Spain 

Sweden 

United Kingdom 

Market 

00 ,i 

0.063“ 

0.045“ 

0.08“ 

0.028 

0.034 c 

0.063 6 

0.022 

0.024 

0.069“ 

0.043 c 

0.032 c 

0.05“ 

01, t 

0.031 

0.083“ 

0.034 

0.011 

-0.011 

0.008 

0.002 

0.006 

-0.003 

0.002 

0 

0.028 

VJi 

0.019“ 

0.015“ 

0.035“ 

0.024“ 

0.024“ 

0.036 6 

0.019“ 

0.018“ 

0.03“ 

0.041“ 

0.019“ 

0.019“ 

01,i 

0.046“ 

0.015 

0.039“ 

0.017 

0 

0.039“ 

0.004 

0.008 

0.028 

0.017 

0 

0 

02, i 

0.079“ 

0.149“ 

0.098“ 

0.103“ 

0.138“ 

0.05 6 

0.131“ 

0.108“ 

0.103“ 

0.13“ 

0.144“ 

0.169“ 

03, * 

0.868“ 

0.886“ 

0.88“ 

0.912“ 

0.909“ 

0.898“ 

0.913“ 

0.924“ 

0.888“ 

0.888“ 

0.91“ 

0.88“ 

Vi 

0.967 

0.912“ 

(MSS 6 

0.879“ 

0.876“ 

0.979 

0.875“ 

0.938 6 

0.891“ 

0.895“ 

0.925“ 

0.824“ 

Vi 

6.989“ 

10.835“ 

10.917“ 

12.093“ 

21.464 6 

8.981“ 

14.539“ 

22.519 1 

12.029“ 

9.602“ 

15.456“ 

13.832“ 


Table A. 4: In sample parameters estimate of the marginal Skew-Student-t AR(1)-GJR-GARCH(1,1) model, 
defined in equation The apexes “a”, “b” and “c”, denote the rejection of the null hypothesis of not significance 
of the corresponding parameter, at different confidence levels 1%, 5% and 10%. 



Austria 

Belgium 

Denmark 

France 

Germany 

Hungary 

Italy 

Netherlands 

Spain 

Sweden 

United Kingdom 

Market 

DGT-ArM 

11.37 

17.35 

15.5 

21.18 

11.68 

21.94 

18.16 

23.5 

13.71 

17.66 

23.68 

20.15 

dgt-ar< 2 > 

24.24 

21.93 

26.15 

12.73 

11.83 

8.88 

11.69 

31.33“ 

10.78 

17.24 

17.25 

13.81 

DGT-ARW 

16.97 

19.88 

22.97 

19.58 

21.01 

22.02 

16.03 

20.27 

18.32 

14.02 

26.52 

16.93 

DGT-ARW 

23.29 

22.37 

29.03“ 

13.18 

11.36 

8.6 

14.99 

31.17“ 

14.12 

18.73 

19.58 

18.09 

DGT-H(20) 

22.42 

16.43 

38.57“ 

12.85 

23.83 

58.7“ 

18.75 

18.87 

16.55 

36.57“ 

36.91“ 

18.3 


Table A.5: In sample Goodness-of-Fit test of 
1%(“), 5%( i> ), and 10%( c ) levels. See also 


Diebold et al. 


11998). Significance is denoted by superscripts at the 


Vlaar and Palm 


(1993) and 


Jondeau and Rockinger 


(2006). 



Austria 

Belgium 

Denmark 

France 

Germany 

Hungary 

Italy 

Netherlands 

Spain 

Sweden 

United Kingdom 

GAScop 

-779.875 

-1518.162 

-1415.213 

-2958.542 

-2861.274 

-634.795 

-2616.335 

-2901.243 

-2218.804 

-2510.746 

-2343.528 

SGASC-L2 

-781.645 

-1523.257 

-1424.369 

-2956.711 

-2869.538 

-622.795 

-2612.144 

-2893.598 

-2224.875 

-2504.244 

-2347.613 

SGASC-L3 

-772.105 

-1508.505 

-1408.366 

-2952.458 

-2853.233 

-616.103 

-2605.672 

-2873.243 

-2216.333 

-2491.122 

-2332.282 

SGASC-L4 

-755.814 

-1495.789 

-1398.815 

-2937.515 

-2839.278 

-597.071 

-2583.14 

-2870.729 

-2198.038 

-2473.934 

-2314.826 

STATcop 

-675.052 

-1394.788 

-1371.802 

-2833.301 

-2774.346 

-614.137 

-2496.057 

-2798.156 

-2081.016 

-2469.562 

-2227.762 

MScop-L2 

-770.958 

-1498.406 

-1415.98 

-2936.215 

-2845.161 

-633.146 

-2595.88 

-2875.751 

-2175.874 

-2500.511 

-2335.623 

MScop-L3 

-774.165 

-1512.3 

-1426.806 

-2952.376 

-2862.936 

-623.068 

-2612.905 

-2895.403 

-2201.949 

-2502.917 

-2341.89 

MScop-L4 

-762.478 

-1503.91 

-1414.824 

-2940.936 

-2857.26 

-611.777 

-2598.457 

-2886.811 

-2199.878 

-2489.773 

-2333.726 


Table A.6: Akaike Information Criterium for the different specification of the dynamic evolution of the dependence 
parameter. For each index, the selected specification is denoted in bold. 



Austria 

Belgium 

Denmark 

France 

Germany 

Hungary 

Italy 

Netherlands 

Spain 

Sweden 

United Kingdom 

GAScop 

0.815 

0.217 

0.134 

0.101 

0.071 

0.42 

0.052 

0.296 

0.275 

0.443 

0.287 

SGASC-12 

0.975 

0.038 

0.13 

0.096 

0.046 

0.419 

0.047 

0.027 

0.122 

0.358 

0.756 

STATcop 

0.987 

0.037 

0.027 

0.139 

0.074 

0.493 

0.083 

0.196 

0.065 

0.039 

0.133 

MScop-L2 

0.935 

0.199 

0.098 

0.099 

0.055 

0.81 

0.027 

0.328 

0.225 

0.048 

0.565 

MScop-L3 

0.923 

0.236 

0.127 

0.013 

0.048 

0.729 

0.309 

0.214 

0.154 

0.027 

0.862 

MScop-L4 

0.902 

0.089 

0.163 

0.011 

0.012 

0.75 

0.103 

0.228 

0.142 

0.013 

0.876 


Table A.7: Test on the correct 
goodness of fit test obtained using 


copula specification. The table reports the p-values of the Anderson- Darling 
1000 bootstrap resamples, see 


Manner and Reznikova 

(2012 

1 and 

Patton 

(2012 
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Austria 

Belgium 

Denmark 

France 

Germany 

Hungary 

Italy 

Netherlands 

Spain 

Sweden 

United Kingdom 

id 

0.081 

0.009 

0.33 

1 

0.001 

1 

0.179 

1 

0.034 

1 

0.005 

a 

0.147 

0.377 

0.712 

0.117 

0.109 

1 

0.657 

0.037 

0.118 

1 

0 

a 

0.909 

0.009 

0.307 

1 

0.005 

1 

0.344 

1 

0.001 

0.995 

0.005 

V 

0.826 

0.721 

0.001 

0.057 

0.303 

1 

0.445 

1 

0.83 

0.681 

0.094 


Table A. 8: Likelihood Ratio p-values for the hypothesis of regime independence of each parameter combination in 
the set {u>, /?, a, uj. Numbers in boxes represent those specifications with the lowest AIC across that with a p-value 
below the significance level of 5% per each country. 



Austria 

Belgium 

Denmark 

France 

Germany 

Hungary 

Italy 

Netherlands 

Spain 

Sweden 

United Kingdom 


0.08459“ 

0.00133“ 

0.00895“ 

0.03873“ 

0.00462“ 

0.0186“ 

0.09411“ 

0.04758“ 

3e - 05“ 

0.17078“ 

0.00059“ 

A 

0.90574“ 

0.99864“ 

0.99385“ 

0.98516“ 

0.99954“ 

0.98316“ 

0.96179“ 

0.98173“ 

0.99998“ 

0.92881“ 

0.99869“ 

ai 

0.02438“ 

0.01183“ 

0.01185“ 

0.04425“ 

0.0017“ 

0.01939“ 

0.0665“ 

0.04423“ 

0.0031“ 

0.0595“ 

0.00975“ 


0.1471“ 

0.08396“ 

0.22526“ 

- 

0.00425“ 

- 

- 

- 

0.05423“ 

- 

0.21305“ 

A 

0.9233“ 

0.98971“ 

0.92492“ 

- 

0.99802“ 

- 

- 

- 

0.99197“ 

- 

0.92514“ 

<*2 

0.07714“ 

0.03378“ 

0“ 

- 

0.00569“ 

- 

- 

- 

0.04082“ 

- 

0.04297“ 


16.55112“ 

40.19022“ 

17.41423“ 

15.584“ 

150“ 

27.13906“ 

15.73629“ 

18.51802“ 

34.39654“ 

11.9017“ 

33.91062“ 

^2 

20.91397“ 

149.99458“ 

149.96147“ 

- 

20.62929“ 

- 

- 

- 

137.29938“ 

- 

8.76538“ 

712 

0.00079“ 

0.15478“ 

0.04524“ 

- 

0.01017“ 

- 

- 

- 

0.1096“ 

- 

0.01788“ 

721 

0.00047“ 

0.02963“ 

0.00639“ 

- 

0.04714“ 

- 



0.04806“ 

- 

0.01838“ 


Table A.9: Parameters estimate of the SGASC-L2 and GAScop models with Student-t copula. The apexes “a”, 
“b” and “c”, denote the rejection of the null hypothesis of not significance of the corresponding parameter, at different 
confidence levels 1%, 5% and 10%. 



Austria 

Belgium 

Denmark 

France 

Germany 

Hungary 

Italy 

Netherlands 

Spain 

Sweden 

United Kingdom 

Pi 

0.42 

0.454 

0.62 

0.861 

0.998 

0.501 

0.841 

0.861 

0.581 

0.832 

0.222 

P2 

0.742 

0.997 

0.903 

- 

0.789 

- 

- 

- 

0.996 

- 

0.888 

Wi 

0.016 

0 

0.052 

0.282 

0.692 

0.005 

0.246 

0.241 

0.004 

0.296 

0 

TO2 

0.086 

0.66 

0.006 


0.125 




0.584 


0.465 

HLi 

7.518 

555.2 

132.628 

69.633 

4783.995 

45.248 

25.924 

56.414 

33868.316 

13.495 

542.058 

hl 2 

11.281 

194.099 

14.419 

- 

477.173 

- 

- 

- 

225.685 

- 

14.034 

Di 

2114.657 

32.75 

155.415 

- 

20.211 

- 

- 

- 

19.808 

- 

53.395 

Dj 

1271.144 

5.461 

21.105 

- 

97.355 

- 

- 

- 

8.124 

- 

54.935 


Table A.10: State identification. For each state l = 1 , 2 ,... ,L, pi, wi denote the long run value of the dependence 
parameter and that of the tail dependence coefficient of the student-t copula, respectively, pi = E | St = l) = 
A ^ j , where A (•) is defined in equation (|9|. HL; is the half life of the GAS dynamics, which can be calculated as 

/ HLA - 1 ral\ n 

the solution of the equation 1 — exp f — /3 l 1 P! j = —^^^ -— 


, D; = represents the expected 


duration (in days) of state l and zui = E (limu-m IP (ms < u \ Uj < u) \ St = l)- 
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Austria 

Belgium 

Denmark 

France 

Germany 

Hungary 

Italy 

Netherlands 

Spain 

Sweden 

United Kingdom 

AE 

1.58 

1.91 

1.65 

1.68 

1.55 

1.2 

1.6 

1.68 

1.29 

1.89 

1.85 

UC 

0.14 

0.04 

0.13 

0.14 

0.18 

0.66 

0.16 

0.12 

0.46 

0.06 

0.07 

CC 

0.34 

0.04 

0.13 

0.32 

0.4 

0.61 

0.15 

0.12 

0.66 

0.17 

0.07 

ADMax 

3.45 

3.51 

3.5 

3.52 

3.52 

3.13 

3.51 

3.52 

3.53 

3.49 

3.52 

ADMean 

0.27 

0.34 

0.29 

0.3 

0.28 

0.19 

0.28 

0.3 

0.23 

0.33 

0.33 


Table A. 11: Actual over expected (A/E) ratios, p-values of the Unconditional Coverage (UC) and Conditional 


Coverage (CC) tests of Kupiec (1995) and Christoffersen (1998), Maximum absolute deviation (ADmax) and Mean 
absolute deviation (ADmean) for the CoVaRf ^ , j — 1,2,..., 11. The AE indicator is calculated as the ratio between 
the realised and expected CoVaR.f j = 1, 2,..., 11 exceedances, while the hitting values of the UC and CC tests 


are calculated using the procedure suggested by Girardi and Ergiin 


(2013). 
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Figure A.l: Marginal empirical PIT distributions. The histogram is divided into 20 bins, red lines represent 5% approximated confidence intervals. For 
more informations seelJondeau and Rockinger (2006) and Diebold et al. (1998]. 
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Figure A.4: Predicted dependence measures: linear correlation (light grey line), Kendall’s tau (grey line), Spearman’s rho (black line), negative tail 
dependence (dotted black line). The black line in the bottom figures represents the predicted marginal volatilities, i.e. &j,t+ i|t, f° r t = 1,2 ,... ,T and 
j = 1,2,..., 11. Vertical dashed lines represent major financial downturns: for a detailed description see Table I A. 1J 
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Figure A.5: Out-of-sample evolution of the ACoVaR (top graph) and ACoES (bottom graph) weekly averages. 
Vertical dashed lines denote the Lehmans failure (September 15, 2008), the peak of the onset of the recent global 
financial crisis (March 9, 2009), the downgrading of the US sovereign by S&P (August 5, 2011), the protraction of the 
Greek austerity package (February 12, 2012), the drop of car sales to the lowest recorded level (September 17, 2013), 
the Italian unemployment rate reaches the record high since the 1977 (November 28, 2014), the Greek government 
unilaterally broke off negotiations with the Eurogroup (June 26, 2015). A complete timeline of the crisis events can 
be found in Table 


A.l 
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